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Abstract 

We  propose  and  analyze  a  remote  atmospheric  lasing  configuration  which  utilizes 
a  combination  of  an  ultrashort  pulse  laser  to  form  a  plasma  filament  (seed 
electrons)  by  tunneling  ionization  and  a  heater  pulse  which  thermalizes  the  seed 
electrons.  Electrons  collisionally  excite  nitrogen  molecules  and  induce  lasing  in 
the  ultraviolet.  The  lasing  gain  is  sufficiently  high  to  reach  saturation  within  the 
length  of  the  plasma  filament.  A  remotely  generated  ultraviolet  source  may  have 
applications  for  standoff  detection  of  biological  and  chemical  agents. 
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The  atmospheric  propagation  of  ultrashort  pulse  lasers  (USPLs)  having  ~  TW  power 
levels  and  pulse  durations  of  ~  100  fsec  are  strongly  affected  by  various  interrelated  processes 
[1,2].  These  include  diffraction,  Kerr  focusing,  group  velocity  dispersion,  spectral  broadening, 
self  phase  modulation,  photo-ionization,  plasma  defocusing,  and  pulse  energy  depletion.  An 
intense  laser  pulse  propagating  in  air  can  be  longitudinally  and  transversely  focused  at  remote 
distances  (>  km)  and  reach  intensities  sufficient  to  ionize  the  air.  Longitudinal  pulse 
compression  is  achieved  by  introducing  a  negative  frequency  chirp  on  the  pulse  while  nonlinear 
transverse  focusing  is  caused  by  the  Kerr  effect  [1-4].  As  a  result  the  laser  intensity  is  enhanced 
(>  10  ‘  W/crrT),  resulting  in  ionization  and  a  plasma  filament  -  1  m  in  length,  ~  100  pm  in 
radius,  and  electron  density  of  Ne  ~  10lf,cm“3  [2]. 
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A  remote  UV  generation  scenario  has  recently  been  proposed  in  which  air  molecules  are 
Raman  or  two-photon  pumped  to  the  appropriate  excited  state  [5,6],  In  this  scheme,  molecular 
excitations  are  achieved  via  simultaneous  action  of  two  synchronized  picosecond  laser  pulses 
with  a  frequency  difference  or  sum  which  is  resonant  with  a  transition.  In  another  publication, 
back  scattered  fluorescence  from  N2  molecules  and  ions  has  been  observed  in  experiments  in 
which  a  USPL- generated  plasma  filament  was  optically  pumped  [7], 

We  propose  and  analyze  a  remote  atmospheric  lasing  configuration  based  on  collisionally 
exciting  the  N2  lasing  line  at  X  =  337  nm  in  a  healed  plasma  filament,  as  illustrated  in  Fig.  1.  A 
secondary  heater  pulse  is  used  to  maintain  and  thermalize  the  electrons.  Thermal  electrons  of  a 
few  eV  can  induce  a  population  inversion  in  N,  leading  to  the  vibrational-electronic  transition 

C’nu  — >  at  337  nm.  Nitrogen  lasers  are  typically  based  on  collisional  excitation  using 

electrical  discharges  [8-10],  In  the  present  configuration  the  UV  radiation  is  generated  in  both 
the  forward  and  backward  directions  along  a  filament  with  a  gain  sufficiently  high  to  reach 
saturation.  For  a  1  p  m  wavelength  heater  beam  the  Kerr  nonlinearity  can  be  the  source  of  third 
harmonic  generation  and  a  seed  for  lasing  at  337  nm. 

The  density  matrix  and  Maxwell  equations  [11,12]  together  with  the  electron  heating  and 
ionization  equations  are  solved.  The  N2  lasing  model  consists  of  an  open  two  level  system 
denoted  by  levels  3  (upper)  and  2  (lower)  respectively.  Since  levels  3  and  2  are  weakly  excited 
the  population  of  level  1  (ground)  is  taken  to  be  fixed.  The  population  of  the  excited  levels  are 


given  by 


3  V, / 3  x  —  vCE _I3  N  +  vct:22  N2  v3  JV3  v,(im  (^3  V2), 


(la) 


(lb) 


where  r  =  t  -  z!  c ,  z  is  the  position  within  the  filament,  t  is  time,  and  c  is  the  speed  of  light. 
Nn ,  are  densities  of  the  «[h  level,  N  ~  =  0.8  Na  is  the  ground  state  density, 

Nu  —  2.7x1 019  cm-3  is  the  air  density,  v2  =  ^CD3I  +  t'ou:  +  +  ^2  is  the  decay  rate  of  level 

3,  v2  =  vCD  2l  +  vC£.ri  +r21  is  the  decay  rate  of  level  2,  IV  is  the  radiative  (spontaneous)  decay 
from  level  i  to  j,  VCF  i}  is  the  collisional  excitation  rate  from  level  i  toj,  and  VCDjJ  is  the 
collisional  de-excitation  rate  from  level  i  to  j  (i,  j  -  I,  2,  3),  vstm  =  axljm  I  Ifico  is  the  stimulated 
emission  rate,  /  is  the  UV  laser  intensity,  co  is  the  lasing  frequency. 
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astim  =  (3  ):  / 4  jt)  (f32  /  /(Agj322  +  y\3)  is  the  stimulated  emission  cross  section. 


ri2  =  (4/3)(o>/c)’  fil2'  I  ft ,  Hn  *s  dipole  moment  between  levels  3  and  2, 

y2i  =  y23l]  +  (1/2)(v2  +  v3)  is  the  dipole  decay  rate,  and  y^is  the  dipole  dephasing  rate  due  to 

collisions  (no  population  transfer).  The  dephasing  rate  y23n  is  the  dominant  contribution  to  y2i , 
The  lasing  intensity  is  given  by 

(a/3z  +  rj)/  =  £/OT,  /(crw.,2)  +  A HI,  (2) 

where  Td  is  the  spatial  damping  rate,  £  is  a  geometric  filling  factor  associated  with  the  seed 
radiation  within  the  plasma  filament,  Istim  —  ftojcN 3,  rfad32  is  the  radiative  lifetime  from  level  3 
to  2  including  the  effects  of  collisions,  and  A  N  =  N3  -  N2. 

For  the  parameters  under  consideration  a  steady  state  population  inversion  is  difficult  to 
achieve.  However,  a  transient  population  inversion  can  be  achieved  if  vCE  13  >  vCF  i2 .  For  a 

constant  heater  pulse  amplitude,  the  population  inversion,  for  N3(0)  =  iV2(0)  =  0,  is  shown  to  be 
AAT(r)  =  N(vCKI3  -  vC£13)(l  —  r/rjt,  (0<r<ri,)  (3) 

where  rm  —  2  ( vt-£|3  —  vC£  |2  )Kvck.  13  (v3  “b  vcD.yi  )  “  vce,n  (v2  vce,23  ))  ’s  typically 

<  1  nsec  and  for  T.  =  2eV  and  Ne  =  10!6cm~3,  the  rates  are  vC£)3  =  2xl05  sec"1  and 
vC£1,  =  1 05  sec-1 . 

The  collisional  excitation  rates  are  dominated  by  electron  excitations  while  the  de¬ 
excitation  rates  are  dominated  by  molecular  collisions.  The  excitation  rate  is 
vCEjj  =  Nc(ijCEi)(v)  vJi  where  aCF:j  denotes  the  excitation  cross  section  from  level  i  to  j  [13]  and 


(  )  denotes  an  average  over  the  electron  Maxwellian  distribution  at  temperature  Tt .  The  life 
times  t,  -  1/v,  and  r2  -  l/v2  are  determined  by  both  spontaneous  emission  and  molecular 
collisions  but  are  dominated  by  molecular  collisions.  Electron  de-excitation  contributions  to  the 
life  times  are  negligible.  In  addition,  since  level  2  is  a  metastable  state,  r32  »  T,, .  The 
lifetimes  of  the  upper  and  lower  levels  are  given  by  [14,  15], 

r3  «  (F32  +  2.2x10*  Pftorr])-1  =  0.6nsec  and  r2  =  (f2l  +  1.9x10s  />[torr])“l  =  7 nsec,  where 


MTn  =  38 nsec  ,  1/f,,  —  9gsec  and  P  is  the  pressure.  The  radiative  decay  time  xraA  32  includes 
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both  spontaneous  and  coilisionally  induced  radiative  transitions  but  is  dominated  by  molecular 
collisions,  i.e.,  zradn  ~  r3.  We  take  yt">l)  to  be  ~10vw,  where  vm  ~  1010sec~‘  is  the  molecular 
collision  frequency  in  air.  Using  the  above  parameters  the  stimulated  emission  cross  section  is 

°sum  ^  10~ucm2. 

The  electron  temperature  is  given  by 

(3/2)— ( =  (J-E„)  -  Qc  Na  Nt{  1  -  TJT) ,  (4a) 

d  r 

(3/2)  Na  =  Q'N'N'O-TJT'),  (4b) 

dr 

where  (j-E,,)  is  the  ohmic  heating,  Qc  is  the  cooling  rate  due  to  all  inelastic  collisions  and  Tv  is 
the  N2  vibrational  temperature.  The  cooiing  rate  Qc  includes  internal  excitation  of  vibrational, 
rotational  and  electronic  states  of  air.  At  low  energies  the  dominant  inelastic  cooling  process  is 
the  excitation  of  vibrational  states  of  N2.  The  cooling  rate  is  obtained  from  the  CHMAIR  code 

[16,17]  and  can  be  modeled  by  Qr  =  3.5x1 0“8  exp(-5/(3r„))  +  6.2 xlO"11  exp(-l  /(3T.))  [18], 
where  Te  <  2eV  ,  Qe  is  in  units  of  eV  cm’/sec ,  and  T  is  in  units  of  eV  . 

The  ohmic  heating  rate  is  =  (ttr  / 8tt) /vf,  where  cop  =  (4 icq2  N^/m)111  is  the 

plasma  frequency,  ve  =  ven  +  vti  is  the  electron  collision  frequency  which  is  the  sum  of  the 
electron-neutral,  vm[sec”']  ~  1CT7  Na  [cm“7]rel/2[eV]  and  electron -ion, 
v^fsec-1]  -  10"5  Ne  [cm'3]r,"3/2[eV]  contributions,  Eeff  =  (ve/0Jg)/(  1  +  v)le%f2  Eu  is  the 

effective  electric  field,  En  is  the  field  amplitude  of  the  heater  beam,  Ia  =  c\Et\~  l %7t  is  the  heater 
intensity,  and  coo  is  the  frequency. 

The  electron  density  is  given  by 

*NJ*T  =  vigaN'-pN,N'-rtNet  (5a) 

a/V^/dr  =  rjNt  -  0_+N_N+,  (5b) 

where  Viim  is  the  collisional  ionization  rate,  /?  is  the  electron-ion  recombination  rate,  /?_+  is  the 
ion  recombination  rate,  jj  is  the  attachment  coefficient,  N_  is  the  negative  ion  density  and 
N+  =  Ne  +  JV_ . 
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The  collisional  ionization  rate  is  vhn  =  vton(N2)  +  vion(02)  where 

v„(X)  =  v,  ( TJUjn(U ,  IT,  +  2)exp(-Ux  IT,)  [19],  X=(AfJ,0!),  [/„,  =15.6  eV, 

UQ t  =12.1  eV,  vN  —  7.6xl011  sec-'  and  =10"  sec-1. 

The  dissociative  recombination  coefficient  [18]  is  /f[cm3  /  sec]  =  1 ,5x10  s  Z’,-0'7  for 
Te  <  0.1  cV  and  /7[cm3  /  sec]  =  2x1 0"R  7y0'56  for  Te  >0.1  eV  .  The  attachment  rate  [18]  is 
^[sec-1]  =  a2  Ng  +  a,  N* ,  where  or2[cm3/sec]  =  2.75xl0"m  Tr~°s  exp (~5/Tf)  and 
£r3[cm6/sec]  =  1.5X10"3'  T~l  exp(- 0.052 fTe).  Typical  values  for  air  [18,20]  at  STP  are 
P  ~  3x10  *  cm3  sec'1,  p_+  ~  2xl0"6cm3sec"1  and  ij  ~  10s  sec-1. 

From  Eq.  (2),  the  lasing  intensity  in  the  transient  regime  is 

I(z,T)  =  e  /f  (exp((rc(T)  -  rd)z  - 1),  (6) 

C  had, 32  (rc;(T)  -  rj 

where  Tc(r)  =  aslim  A/V(r)  is  the  spatial  growth  rate  and  AV(r)  is  given  by  Eq.  (3). 

Atmospheric  lasing  is  simulated  by  solving  Eqs.  f  I ),  (2),  (4),  and  (5)  using  a  heating 
pulse  with  intensity  profile  lluaKr{z,t)  =  1o  sin:(jrr/r0),  0^ r5r(, ,  and  wavelength  Xn  =  1pm  . 
Ionization  by  the  USPL  is  assumed  to  have  created  a  filament  with  a  uniform  density, 

Ne  (0)=1016  cm-3  at  T.  (0)= 0.025  eV.  We  hold  the  fluence  of  the  heater  pulse  constant,  i.e., 

I0T0  =  285  J/cnr ,  vary  the  pulse  duration,  and  observe  the  growth  of  the  UV  radiation  along  the 
filament.  The  geometric  factor  is  e  ~  (rf  !2La)2 ,  where  rf  ~  100pm  is  the  filament  radius, 

Lg  =  1/FC  is  the  gain  length,  and  Trf  =  0 . 

Figure  2{a)  plots  the  UV  fluence  as  a  function  of  position  along  the  filament  for  heater 
pulses  with  xtt  =  0.3,  0.4,  and  0.5  nsec  .  Initially  there  is  exponential  growth,  which  is  followed 

by  linear  growth  as  the  population  inversion  is  modified  by  the  UV  radiation.  Figure  2(b)  plots 
the  gain  length  of  the  fluence  as  a  function  of  position.  The  gain  length  during  the  exponential 
gain  regime  is  dependent  on  the  heater  pulse  intensity,  i.e.,  Lc  =  0.5,1,  and  1.5  mm  for 

ro  =  0.3,  0.4,  and  0.5  nsec  ,  respectively.  The  gain  lengths  are  in  good  agreement  with  those 
used  in  Eq.  (6).  Exponential  gain  ends  after  ~13Z,C  in  all  three  cases. 
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Figure  3  plots  the  temporal  profile  of  the  UV  radiation  and  population  inversion  at 
various  positions  along  the  filament.  In  Fig.  (3),  the  heating  pulse  duration  is  to=  0.4  nsec ,  the 

heater  intensity  is  In  =  7x10"  W/cnr ,  and  the  duration  of  the  initial  AN  is  in  good  agreement 

with  Eq.  (3).  Figure  3  shows  a  -  0.1  nsec  UV  pulse  forming  near  the  back  of  the  inversion 
region.  As  the  UV  pulse  grows  AN  is  depleted  and  the  peak  of  the  UV  pulse  moves  forward 
relative  to  the  heater  pulse.  After  ~  2  cm  of  propagation,  the  peak  UV  intensity  is  ~  8  MW/cm2 
and  the  fluence  is  -  0.5  mJ/cnr. 

We  have  analyzed  a  lasing  scheme  that  can  produce  UV  in  the  atmosphere  at  kilometer 
distances.  Other  molecular  transitions  may  also  be  excited  and  result  in  lasing,  e.g.,  N2  and  02 . 
Non-uniformity  of  the  plasma  filament  may  also  affect  the  gain  and  saturation  length  [21]. 

Acknowledgements.  The  authors  are  grateful  to  Drs.  R.  Femsler,  S.  Suckewer,  and  A.  Zigler  for 
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Figure  1 :  Schematic  diagram  of  remote  lasing  configuration.  An  ultrashort  pulse 
laser  creates  a  plasma  filament  of  seed  electrons  that  is  heated  by  a  secondary 
heater  pulse.  In  the  energy  level  diagram  for  N2  (right),  the  energetic  electrons 
collisionally  excite  the  N2  molecules  and  induce  lasing. 
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Figure  2:  (a)  UV  fluence  as  a  function  of  propagation  distance,  z,  for  heating 
pulses  with  ka  =  1  fim,  lnxn  -  285  J/cm2  and  pulse  durations 
r„  =  0.3,  0.4,  and  0.5  nsec .  (b)  Fluence  gain  length,  Lq,  as  a  function  of  z, 
corresponding  to  panel  (a). 
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Figure  3:  UV  intensity,  /  (solid  curves)  and  normalized  population  inversion 
AN  /  N  =  (N3  —  N2)/  N  (dashed  curves)  as  functions  of  time  at  (a)  z  —  0.9cm, 

(b)  z  =  1.4cm,  and  (c)  z  =  2cm,  for  a  heating  pulse  with  X0  =  1  fmt, 
l0z0  =  285  J/crrf  and  pulse  duration  rn  =  0.4nsec 
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